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Abstract 

In this paper, we present a concept of the averaged probability density function (APDF) for studying 
compressible turbulent reacting flows. The APDF is defined as an ensemble average of the fine grained 
probability density function (FG-PDF) with a mass density weighting. It can be used to exactly deduce 
the mass density weighted, ensemble averaged turbulent mean variables. The transport equation for 
APDF can be derived in two ways. One is the traditional way that starts from the transport equation of 
FG-PDF, in which the compressible Navier-Stokes equations are embedded. The resulting transport 
equation of APDF is then in a traditional form that contains conditional means of all terms from the right 
hand side of the Navier-Stokes equations except for the chemical reaction term. These conditional means 
are new unknown quantities that need to be modeled. Another way of deriving the transport equation of 
APDF is to start directly from the ensemble averaged Navier-Stokes equations. The resulting transport 
equation of APDF derived from this approach appears in a closed form without any need for additional 
modeling. The methodology of ensemble averaging presented in this paper can be extended to other 
averaging procedures: for example, the Reynolds time averaging for statistically steady flow and the 
Reynolds spatial averaging for statistically homogeneous flow. It can also be extended to a time or spatial 
filtering procedure to construct the fdtered density function (FDF) for the large eddy simulation (LES) of 
compressible turbulent reacting flows. 

1.0 Introduction 

The advantage of PDF method in the turbulent reacting flow simulations has been studied and 
examined (Refs. 1, 2, 3, and 4). Here we present the concept of an ensemble averaged probability density 
function (APDF) that may be considered as an extension of the conventional probability density function 
(PDF) method for dealing with compressible turbulent reacting flows. We do not assume that the mass 
density is a function of temperature and composition only, because this is only valid for the Mach number 
being near to zero (see Ref. 4). The APDF is defined by an ensemble average of a mass density weighted 
fine grained probability density function. It can be used to exactly define the mass density weighted 
ensemble averaged turbulent mean variables, very similar to the way in which the conventional PDF 
defines the turbulent mean variables. 

Two ways of deriving APDF transport equation will be described in detail. One is to start from the 
transport equation of the fine grained probability density function (FG-PDF), in which the compressible 
Navier-Stokes equations are embedded. By taking ensemble average over the equation of FG-PDF, the 
resulting transport equation of APDF in its traditional form contains conditional means of all terms from 
the right hand side of the Navier-Stokes equations, except for the chemical reaction term. These 
conditional means are new u nk nown quantities that need to be modeled. Another way of deriving the 
transport equation of APDF is to start directly from the ensemble averaged compressible Navier-Stokes 
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equations. Using the relationship between APDF and the mass density weighted, ensemble averaged 
turbulent mean quantities, and under a “conservational” condition, the derived transport equation of 
APDF results in a closed form without any need for additional modeling. No terms of conditional means 
will appear in this conservational form of APDF equation, except when a marginal APDF equation such 
as the scalar APDF is considered. 

The methodology of ensemble averaging procedure presented in this paper can be extended to other 
averaging procedures: for example, the Reynolds time average for statistically steady flows and the 
Reynolds spatial average for statistically homogeneous flows. In these cases, the Reynolds time (or spatial) 
averaged probability density function will be defined, which then defines the mass density weighted 
Reynolds time (or spatial) averaged turbulent mean variables. 

It should be noted that for the constant density turbulent flow, the APDF will reduce exactly to the 
conventional PDF if the expectation of FG-PDF can be viewed as the ensemble average of FG-PDF. The 
conservational PDF transport equation has been reported in Reference 5. 

The methodology described in this paper can also be extended to a time or spatial filtering to 
construct the transport equation for a mass density weighted filtered FG-PDF for large eddy simulations 
(LES). In this case, the mass density weighted filtered FG-PDF is named as “FMDF” or “DW-FDF” (see 
Refs. 3 and 6), which is a random quantity and can be used to define the mass density weighted filtered 
turbulent variables to describe large scales of compressible turbulent reacting flow. 

Finally, it is noted that the equations for various (PDF, APDF, DW-FDF) methods are sought to be 
similar, so that the solution procedures would be similar regardless of the solution variables. 

2.0 Ensemble APDF 

In this section, we will use the FG-PDF to define the ensemble APDF, and then explore the 
relationship between APDF and ensemble averaged turbulent mean variables. This will provide the basis 
for establishing the transport equation of APDF. 

2.1 Fine Grained Probability Density Function fh (V ; x, t), /$ (i|/; x, t) 

The FG-PDF for turbulent velocity and scalars (e.g., species mass fraction or composition, internal 
energy) are defined as follows (Pope (Ref. 1 )), 

3 

ti{V-x,t) = h(U{x,t)-V) = Y\K U i (*> O-K) ( 1 ) 

1=1 

M+l 

/i (l|/; AT, 0 3 8 (<£ ( x, 0 - V|/ ) 3 Y[ 8 (° m (*> 0 - Vm ) (2) 

m = 1 


where 8 denotes the delta function, U(x,t) is the turbulent (random) velocity vector (U u U 2 , U 3 ), 0(.v,t) is 
the turbulent (random) scalar array ( , , O 2 , • • • , O M , ® M ( l ] , for example, M species mass fractions and 

one internal energy © M+1 = e , the x,t denote the physical space variable (x t , x 2 , x 3 ) and the time t, V= 
(V h V 2 , V 3 ) and i|/ = (v)/, , i|/ 2 ,- • ■ , Va/+i) arc the sample space variables for IZ(xJ) and 0(x,t), 

respectively. 

2.2 Ensemble Averaged Turbulent Variables and APDF F U (V; x,t ), F 0 (\\i; x,t) 

We start with the definition of ensemble averaged turbulent variables and the definition of ensemble 
averaged probability density function APDF. And then explore the relationship between the two. 
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2.2.1 Definition of Ensemble Averaged Turbulent Variables (j)(.v,t) , <|) (. x,t ) 

In the case of compressible turbulent reacting flow, we often deal with two types of ensemble 
averaging: one with the mass density weighting, the other without the mass density weighting. The 
ensemble averaged turbulent variable without the density weighting is denoted by <|)(jc,t) and is defined 
as 


<k*,o= iim 


(3) 


n = 1 


where (j) <n) is the n-th independent realization of turbulent variables, for example, the velocity 


components u\ n) , the dcnsityp'" ) , pressure P' n> , the species mass fraction and the internal energy 


,(«) 


3<») 


(«) 


M 


e (n> — ^ . N is the total number of independent realizations. The mass density weighted 


m = 1 


ensemble averaged turbulent variable is denoted by ()> (x, t ) and is defined as 


<K*>0 = ^ 


(4) 


These ensemble averaged variables 4 >(jc,?) , §(x,t) represent the turbulent mean variables and are 
measurable in experiments. They only contain relatively low frequency and low wave number 
parts of the turbulent motion comparing with the non-averaged turbulent variables <t>(x,t). 

2.2.2 Definition of Ensemble APDF 

We define the following ensemble averaged probability density functions APDF: 


F,j(V- x,t) ^ Urn -^£p (n) (x,t)-V) 

1 


n = 1 
N 


(5) 


F 0 ( W,x,t)= lim — V p ( ' !) (.v,t)8(o ( ' !) (.v,t)-v|/) 

jV— >00 Tv V ' 

n = 1 


Obviously, thus defined F u (V; x,t) or F (l , ( \|/; x,t) is no longer a random quantity. And they satisfy the 
following “normalization” property: 


00 00 . N N 

| F v (V; x,t)dV = | lim (x,t)h(lJ^ (x,t)-V}dV = lim — ^ p (n) (x,t) - p 


iv-> oo N , 

-00 n = 1 

N 


n = 1 
N 


( 6 ) 


~ ~ iV 

| x,t)dif / = | lim — ^P (H> (v,t)8^0 (,i> (v,?)-v|/j<f^= lim — ^p (n) (.v,t) = p 

-oo -oo ^°° »= 1 N->CO n= j 


Where the integration proceeds in the entire sample space of 17 (- 00 , 00 ) and i)/, (0,1) . 

2.2.3 Relationship Between APDF and Ensemble Averaged Turbulent Variables 

With the definition of APDF described in Equation(5), we can exactly deduce the ensemble averaged 
turbulent variables that are defined in Equations (3) and (4). For example, 
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£Vf„(F; x,t)dV - \y( lim (x, >)s(t/« (.v,,)- r) 


AT-*» jV ' . 

V n — 1 




( i ^ 
«=i 


W 


n = 1 

= pi/ = p(x, ?)£/(*,?) 


where the sifting property of FG-PDF has been used in the second line of Equation (7). 
Similarly, we have 


(7) 


J vF 0 (y; x,t)d\\i = p® = p(x,t)®(x,t) (8) 

If we denote the left hand side of Equations (7) and (8) as a “mean” operation using APDF: (U^ and (<J>) , 
we may write Equations (7) and (8) as 


(u) = J += ° VF v (V ; x,t)dV = pU = p(x,t)U(x,t) 

(9) 

(<J>) = J ^ v|/P (t (\|/;.v,t)^-p(D = p(.v,t)a)(A:,0 

(10) 


For any function Q{U {x,t)) or W (0(x,t)) , it is easy to verify that 


{Q(U)) = ^Q(V) F u{V^xS)dV = p(x,t)Q(U{x,t)) 


( 11 ) 


(w (<£)) = J ^ W (i| /)F 0 (\|/; x,t)dy = p (x,t)W(<X>(x,t)) 

Furthermore, we may consider the derivatives VP, Vfj',VO as new random quantities and legitimately 
write 


(VP) = p VP, (VU) = pVU, (VO) = pV<D (12) 

However, because of the variable density, the mean ) does not have the differential commute property, 
i.e., 

(vp)^v(p), (VU)*V(U), (vo)^v(o) ( 13 ) 

because p VP ^v(pP), p V£/^ v(pi/), pVO^v(pO). 

It can be verified that the following “second moment” is true 

(u,Uj) = \ + y< Vj F v (V ; x,t)dV = p U t Uj = p(x,t)Up' j {x,t) (14) 
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(15) 


We can also write for the joint variables as 

(UjOj ) = | Vj \|/,- F u <t) (I V , v| /; *, 0 c?FA|/ = p (x,t ) D,® , (x, t) 
where, F v 0 (F, v|/; x,t) is the velocity-scalar joint APDF defined as 

F u ^ v ’V’ x ’ t )= Km ^X p(H, (- v »0 8 ( t/( " > ( A: ’0 _f/ ) 8 ( o( ” ) (- v ’0 _v i / ) ( 16 ) 

n = 1 

Here, the velocity-scalar joint fine grained probability density function / 6 ', (|) is introduced: 

/^®(F,v;jc,O = 8(^(^0-F)S(<I>(^0-v) (17) 

The joint APDF ^(F,^; x,t) also satisfies the normalization property: 

00 00 N 

f [ 7v,a,(F,i|/; x,t)dVd\y= lim — V p ( " } (x,t) = p (18) 

j j ’ n — >oo tv 

—00 —00 n = ' 

2.2.4 Conditional APDF and Its Conditional Mean 

From Equation(16), we may follow Pope (Ref. 1) to define a “conditional” APDF on the condition 

® = i|/ as 


F u\t><yW>x, t) 


F u,<s>(yw x,t) 

F^',x,t) 


(19) 


and it has the following property, 


oo 

{ F u\<t> (F| v|/; x,t)dV = 1 

-oo 


( 20 ) 


The “conditional mean” will be defined as 

/ \ r +co 1 r +o ° 

(U(x,t) v|/}= [ VF m (V v|/; x,t)dV = — -[ V F v <j,(F,y|/; x,t)dV 

' ' 1 F. ( ur v / w-c o 


F< 5 (v; x,t) j - a 


F lim — ^p ( '9 (x,t)S(u t ' n) (x,t) - V )§(® ( ' l) (x,t) - v|/ ) 

- d> _a> l^"* 00 A „ =1 

Then we have 

lim — ^p (n) (.*,?)I7 ( " ) S(® ( ' ,) (.*,?) - v|/) = J ^ V F v j0 (F,\|/; x,t)dV = ■ (U(x,t)\ v|/ 


w 


TV— >co TV 


(21) 


( 22 ) 
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And the total mean should be defined as 


Loo 

1 N 

= J + '"J +C °F lim — ^p ( ' ,) (.v,t)8(t/ ( ' i) (A:,t)-F)8(o ( '' ) (.v,t)-v|/)^Fi/\|/ (23) 

1 N __ 

= lim — ^p (n) ( x,t)U (n) (x,t) = p U - p[x,t)U(x,t) 

Equation (22) can be extended to any other turbulent quantities, for example, VP , VU , VO , 5) (O) : 
lim — ^p (n> (V j P) ( " ) 8(f ( " ) (*,t) -V) = F (J ■ (VP\v) 

lim — ^p (n> (Vt/) ( " ) 8(t/ ( " ) (x,t) -V) = F U - (VU\V) 

"n (24) 

Km -|-]r p ( "> (V0) ( ' !) 8(0 ( ”» (x,t) - v|/) = F* • (vo|v|/> 

n = 1 
n = 1 


where VP, VU, VO are viewed as new random variables in addition to P , U, and O . 

3.0 Traditional Form of APDF Equations 

In this section, we perform the traditional way of deriving the transport equation for APDF, i.e., we 
start from the equation of FG-PDF, and the resulting equation for APDF contains conditional means that 
require closure modeling. 

3.1 Traditional Form of F v (V; x,t) Equation 

The transport equation for the fine grained PDF with the variable density p(.v,t) can be written as (see 
Ref. 6): 


dp fu | dp U jfu 
dt dxj 



f 


dP + d 
v dx t dxj v 


2p5, 




ij°kk 


(25) 


where the compressible continuity and momentum equations (31), (32) have been applied. Applying the 
ensemble average operation. Equation (5), on Equation (25) and using the sifting property of FG-PDF, we 
may obtain 
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8F, 


dt 


dFr 

] dx 


1 N 

lim — Y p ( ") 

dV- I 


u+v r u 


n = 1 


^1 gp^ ( " ) 
v p y 


U) 


(«) 


a 

5K. 


lim — Vp ( " 
N^N^i 

n = 1 


JV 


/ , N \(») 
'ij/ 2 v 


(fu) 


(") 


According to Equation (24), the APDF equation for the velocity can be written as 


d JlL +V . d IlL 

dt 7 dX : 


dV; 


\-F 71^ 
b \P Sxi 


V +Fu-(~ 


1 d 


2 ") 
, jj — ~q8 ijS/ik 

P dx A 3 j 


( 26 ) 


At this point, the velocity APDF equation is general but unclosed because of the u nk nown terms of the 
conditional means. The last term in Equation (26) is referred as the molecular mixing term in which 
significant modeling efforts have been devoted over several decades, nevertheless, major issues still 
remain. Later, we will show that when using a different way to derive the APDF equation, the molecular 
mixing term is in a closed form, and the pressure gradient term is also essentially closed. In addition, we 
noticed that the molecular mixing term is of the order of 0(1/Re), and may become vanishingly small 
when comparing with other terms in the equation in the case of large Reynolds number (Re). 


3.2 Traditional Form of F^iy; x,t) Equation 

Applying the same procedure (in Section 3.1) on the following FG-PDF equation for scalars, 


dp fo , d P u jf<t> _ d , 

' + _ 1 J<s> 


f d r 


\ dx j 


dt dxj d\\i t 

we obtain the APDF equation for the turbulent scalars: 


pr''>^i 

8X J J 


+ p S f (©(x, 0) 


(27) 


dF<& 

dt 


dx,- 


■— P r(0 — 
di|/ - \p dx - dx • 


\v) + F (J) -S i (^)\ i = + 1 (28) 


where S M+l (v|/) = 0 . Equation (28) is also unclosed because the conditional means are unknown. Later, 

we will also see that the molecular mixing term is closed when using a different way of derivation. Again, 
it is of the order of 0(1/Re) and could become vanishingly small at large Reynolds numbers. 


3.3 Traditional Form of Joint APDF F v <J) (V,\\r, x,t) Equation 

We may obtain the velocity-scalar joint APDF equation from the joint fine grained probability density 
function (|) defined in Equation (17). Similar to the relationship Equation (25), we may write 
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dpUjfu^ _ fi I , 
dt dx, dvA U,<5> 


~ dv \ fu,o 


r dp u, dpu^ij ' 


dt 


dx 


J J 




\fi 


U,<5> 


f dp<$> k + dpu j (i) k ' 


dP 


dt 


dx 


J ) 


d^i 


f’ 
i Ju 


^ fix,- dxj 


'fi ' 


2pS u --pS,,A 




c/,® 


K dx J 


pr 


■(*•) 


50, 

j 


l] 3 

+ p5' jt (O(x,0) 


Applying the ensemble average operation (5) on Equation (29) we obtain 


^U&_ + y 


dt 




d 

dV t 

\ F 

\ 1 T 


~Ft 


1 dP 


u,<b 


O dx, 


V,W)+F 1 


1 5 


c/,o 


1 5 r a> dq>k 
U,(S> ^ p fix- 5x - 


P dx y - v 

F,\|/ ) + F U (S> ■ S k (v|/) 


2p5, --pS,S, 


3 


ij^kk 


v,v 


(29) 


(30) 


It is easy to verify that Equation (30) will reduce to Equations (26) and (28) by an integration over the 
sample space i|/ and V, respectively. The equation (30) is general and is not restricted to low Mach 
numbers. 


4.0 Conservational Form of APDF Equations 

In this section, we first review the ensemble averaged compressible Navier-Stokes equations, which 
are unclosed due to the unknown second order moments. Then we use the relationship between APDF 
and the first and second order moments described in Section 2.2 to express the ensemble averaged 
compressible Navier-Stokes equations in terms of an integral equation of APDF. Under a sufficient 
“conservation” condition, the APDF transport equations are established, which are in a closed form 
without additional modeling, unless a marginal APDF equation is considered. 


4.1 Ensemble Averaged Compressible Navier-Stokes Equations 

The Navier-Stokes equations for compressible reacting flows can be written as 

a P ! a P U J Q 

dt dxj 


dp U i d p UfUj 
dt dxj 


dP_ d f 
fix,- dxj V 


2pv(*SE 


--S yS kk ) 


dpe dpU t e 
dt fix,- 


+ PS kk + 2pv 
fix,- 


S,S H — s, s, 


v v 


iPkk 


+ Q 


dpQ„, | fipU,Q„, 
dt fix,- 


fix,- 


pr 


(m) 


50 „ 


/ v 


fix, 


+ W m 


i J 


m - 1,2 


(31) 


(32) 

(33) 

(34) 
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p = P r Yj~ 


m = 1 




c t, = -p k <a 


ST 

dX: 


M 


2 > 

m=l 


(m) 


3X: 


( 36 ) 


Applying the ensemble averaging on the above equations, we obtain the following Reynolds averaged 
Navier-Stokes equations: 


dp , dp Uj _ Q 

dxj 

dpUj dpUjfj _ 8P_ d 
dt dxj dx i dxj 


2pv(5,-is^) 

V 6 J 


dpe dpU t e 
dt dXj 


dd ■ 

~~~ + PSkk + 2p v 

OX: 


SySy —S n S kk 


+ Q 


0p® m dpU i Q> m 
dt dxj 


_d_ 

dX; 


pr 


(m) 


dO„ 


dX: 


+ Wn 


i J 


m = 1,2 


where 


M 


P = P R Y, 


o.„r 




or 

P R y 1 


( 37 ) 

( 38 ) 

( 39 ) 

( 40 ) 

( 41 ) 


9i 


dT 

dX; 


Zpr 1 ”'*, 


m = 1 


dO, 

dX: 


( 42 ) 


In the above equations, k, v and T im) are the molecular heat conductivity, kinematic viscosity and the m-th 
species diffusivity, they have the same dimension (i.e., velocity • length). It is also commonly assumed 
that r ( "° is the same for all species ®,„. The h m , T are the enthalpy of m-th species and the temperature, Q 
is the radiation rate, W m = pS m is the chemical generation rate of m-th species, ® M+1 represents the 
internal energy e, R is the universal gas constant. These equations are general; however, unlike the 
constant density flows, further approximations for the terms on their right hand side are required in order 
to complete the ensemble averaging process. One of such approximations leads to 


2pv 


Sy --5 yS kk 


dUi dUj 
+ - J 


\ dx J 


dx. 



dU k ' 

dx k , 


( 43 ) 
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In which, we have basically neglected the variations of p and p during the averaging process, the value of 
p will be considered as the function of P,T,---. Alternatively, 


Similarly, 


2pv 




'dp Ut , dpUj 

dxj dx i 


2 5p U k " 

3 ij dx k j 


dT 

dX: 


M 




m = 1 


50 

>u 

dX: 




M 


50 „ 


dx, 


m = 1 


8X; 


(44) 


(45) 


p(m) dQ,„ „ p(m) dpQ;» 
dXj dXj 

Furthermore, invoking the turbulent kinetic energy dissipation rate: 


2pv 


SySy-~S u S kk 


= ps 


Where v, k , c v , and F 1 '" 1 are now considered as functions of P,T,--- . 

Then, the ensemble averaged Navier-Stokes equations can approximately be written as 

dp | dp Uj _ Q 
dt dx j 


dpUi 5 p UjU : dp 5 

- + - = + ■ 


dt 


dX: 


dXj dxj 


^dpUi + dpUj _ 2^ dpUk ^ 
^ dx j dx t 3 iJ dx k 


dt dX: dX: 


5pO„, 5p(7,O m _ 5 


dt 


dX: dX: 


\ 


r (m) 5pO 


V 

M 


dx,- 


+ p S m m = 1, 2, • • • , M 


i j 


OJ P^vO,/ 






m = 1 C v m = 1 


or 


- = pi? y^0„,0 M+ i 


q, = ~ 


dpT 

dX: 


M " 

Ipr^ 


i m = 1 


dX: 


(46) 

(47) 

(48) 

(49) 

(50) 

(51) 

(52) 

(53) 
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These equations are still considered as quite general, because i) they are exact if the flow becomes 
incompressible, ii) all the approximations made in Equations (49), (50) and (51) are related only to the 
molecular diffusion terms that are less important and even become negligibly small (1/Re) comparing 
with the convection terms on the left hand side for turbulent flows at large Reynolds numbers (see 
Tennekes & Lumley (Ref. 7) and Pope (Ref. 1)). For the conventional unsteady Reynolds averaged 
Navier-Stoke simulation (URANS), Equation (48) to (51) are often used together with the further 
approximations for Equations (52) and (53): 


M 




<D T 


m = 1 W m V 


pRT 

M 


p RT 1 _ y <E„, y 


q, = - K u 


dpT 

8X; 


(54) 

(55) 


The momentum flux pUJJ . , the energy flux p Up and the species flux p U <X> m are considered to be 

critically important in URANS simulations and should be carefully modeled. Many models in the 
literature, from the simplest zero-equation model (Ref. 8) to the complex two-equation models (Refs. 9 
and 10) have been suggested. In the following Sections 4.2 and 4.3, we will derive the APDF equations 
directly from Equations (48) to (53). It is important to note that in the joint PDF method, the terms 
corresponding to these important turbulent fluxes are always in closed forms with no need of modeling. 


4.2 Conservational form of F u (V;x,t) equation 

Using Equations (9) and (14), the left hand side of Equation (49) can be written as 


dp Ui dp UjUj _ d 


d 


dF rr „ 8F n 


= - I V i F u dV + — f V i V,F u dV= f VJ ^L +Vj — U 
8t dx: dt J Ox, J l 8t dx. . 

J —OO J —00 —00 \ J / 


dV 


(56) 


The pressure gradient term can be written as (via Eq. (52)) 


dP __ d 
dx t dXj 
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V v m = 1 m J 
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m = 1 C v W m 


dV 


00 


V- 


dV k dx k 


M 

z- 


R 


1 V)-Fu 


dV 


(57) 


During the above arrangement, no approximations have been made other than the integration by parts and 
a zero integration property like Equation (59). Note that the pressure gradient term can be in a closed 
form only when the joint APDF, F v ^ , is considered, otherwise, a conditional mean will be unavoidable. 
Now the molecular mixing term in Equation (49) can be written as 
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and further rearranged as 
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J dx, 
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dV- \v t -*~ 


r d 2 V j F u 
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dy+ i r ‘h>£; 


r d 2 V m Fu 

V 8V k dx m J 


dV 


Note, on the third line of Equation (58), we have applied the integration by parts and the following type of 
zero integration similar to the one used in PDF formulations (Pope (Ref. 1)): 


°° d 

\ w m Fv )dv= o 

-oo ^ 


(59) 


Finally, collecting the terms that are in the integrands and factored by V t , we obtain the following 
transport equation for F v (V;x,t) by setting the sum of the integrands to zero, which will sufficiently 
satisfy the Navier-Stoke equations (48), (49) and (52). We refer to this condition as the “conservation 
condition”. 


8 F u. + v 8F u 
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(60) 
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In Equation (60), all the terms are closed, except for the one involving the conditional mean originating 
from the pressure gradient. It should be noted that the unclosed term in the momentum equation, 

i.e., pUjU : , now is closed in the corresponding APDF equation while the less important molecular 

diffusion terms remain in closed forms. 


4.3 Conservational Form of (v|/; x, t) Equation 

Similarly, we may obtain the APDF equation for the scalars (i.e., species and internal energy) from 
Equations (50), (51) as follows: first, we write the terms on the left hand side of Equation (51) as 


5pO,„ 8 f ^ , f dF* J 

-ir = aJ^ F ^ = J^-ar^ 


(61) 
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(62) 


Then, the terms on the right hand side of (5 1 ) can be written as 
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S(S t F Q ) 


d\\i 


(64) 


Where in Equations (63) and (64), we have applied the integration by parts and zero integrations similar 
to Equation (59). Collecting all the integrand terms that factored by \\i m , and invoking the sufficient, 
conservation condition, we obtain 
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(65) 
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This equation also represents the equation of internal energy when S M+l (i|/) = 0 and other source terms 
in Equation (50) are neglected. Equation (65) is essentially closed if we consider the joint APDF F v <t) . 

For the marginal , the convection term is not closed because of the conditional mean . Then, 

this critically important term, corresponding to p C/,-0 in Equation (51), must be carefully modeled 
while the less important molecular diffusion term remains in a closed form. In addition, we noticed that 
the equally important chemistry source term p5 m (0) in Equation (51), which involves complex 
processes of turbulence-chemistry interaction and is very difficult to be modeled accurately, now is closed 
in the APDF equation with no need of modeling. This direct calculation of turbulence-chemistry 
interaction is one of the unique features of PDF methodology. 

4.3.1 Modeling of Unclosed Terms 

Consider the approximation described by Equation (54), i.e., P « pR T /m , we may write from 
Equation (57) 


dP ~ 
dx, dx. 
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=<® M+ i\ V)-Fi 


u 


dV 


( 66 ) 


For the marginal APDF equation of F v , we must model the u nk nown “conditional” mean ^® M+1 . 

While seeking its physical based model the roughest one could be suggested as follows: 

(® M+X \V)*® M+ X=e (67) 

As for the term containing conditional mean, F (X) ■ (U t ■ | \|/^> , in Equation (65), we may start from a more 
general model for the term pUf P m (see Refs. 1 1 and 12): 


p[/ / o w — pU i&m r 


(m) 5p®« 


dx, 


-r 


(m) 


-I Ci 
8 


1 Sij + C 2 Qij 


, op®„ 

Gxj 


( 68 ) 


Where c, , c 2 are model coefficients and were temporally set to c, = c 2 = - 0.24, r 1 /" 1 is the turbulent 
diffusivity for species and is assumed to be the same for all species. This will lead to the following 
model by directly applying Equations (8) and (22) described in Sections 2.2.2 and 2.2.4: 
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( 69 ) 



4.3.2 Summary 

With the models given by Equations (66) and (69), the marginal APDF equations for F v (V;x,t) and 
F* (v|/; x,t) can be written as 


dFjj dF v d 2 I” F Om+i „ ] d f dFjj 

dt 1 dxj dV k dx k [ c v M J dxj ^ dxj 

d f d d (l d dF v 

dV k y dXj dx k J dV k \2> d Xj dx m 

or (70) 

dF u lV dF u_ d 2 R®m + i F d ( d ^dF u 

dt 7 dxj dV k dx k c v M u dV k ^ k dxj dxj 

d f d d (l d dF v 

dV k y dxj dx k J dV k [ 3 dxj dx m 


and 



It can be verified that the APDF Equations (70) and (71) can exactly deduce the ensemble averaged 
compressible Navier-Stoke equations (48) to (51) with the assumed approximations of (66) and (68). 
Elowever, the models described by Equations (66) and (69) are by no means unique. Furthermore, the 

variables ( O 1 /+ 1 , U , , Sy and Q ;/ ) are considered as available quantities during the solution procedure of 
APDF equations. 
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4.4 Conservational Form of Joint APDF F U (il {V ,vf,x,t) Equation 

The transport equation for the velocity-scalar joint APDF F U (b ( V , v|/; x, t ) can be constructed from 
Equations (60) and (65) as 
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It is easy to verify that Equation (72) or (73) will reduce to Equations (70) and (71) by an integration over 
the entire space vj/ and V, respectively. It is noted that the transport equation for joint APDF 

/•CO 

F ( j 0 (V, O; x, t ) is in a closed form; however it contains an integration of ^m+i^u ® ^ rom t ^ ie 
pressure gradient term. 


5.0 Concluding Remarks 

We have derived the traditional form of APDF equations by starting from the equation of FG-PDF. 
The resulting equations contain the “conditional means” for all the terms on the right hand side of the 
embedded Navier-Stokes equations except for the reaction term (see Equations (26) and (28)), which need 
empirical models to make the equations closed. 

We have also derived the conservational form of APDF equations directly from the ensemble 
averaged compressible Navier-Stokes equations using the relationship between APDF and the ensemble 
averaged turbulent variables. In this way, the resulting joint APDF equation is in a closed form, see 
Equation (72) or (73). 

APDF equations both in the traditional and in the conservational forms derived in this study do not 
subject to the low Mach number limitation. 

In the marginal APDF equations, (60) and (65), there are two unclosed terms, one is related to the 
pressure gradient and the other is related to the scalar convection. They are generally much more 
important than the molecular “mixing” terms, especially at large Reynolds numbers. The models 
proposed in Equations (66) and (69) should be further developed and evaluated. 
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It should be pointed out that the methodology used in the present study can be extended to other 
averaging operations: for example, the Reynolds time (or spatial) averaging. Furthermore, it can also be 
extended to a filtering operation to derive the filtered density function (FDF) for large eddy simulation 
(LES) such as the mass density weighted filtered density function DW-FDF for LES of compressible 
turbulent reacting flows (see Ref. 6). 
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